A&A manuscript no. 

(will be inserted by hand later) 

Your thesaurus codes are: 
02.08.1;02.13.2; 11.11.1; 11.19.2 

Global three-dimensional simulations of magnetic field 
evolution in a galactic disk 

II. Gas rich galaxies 

S. von Linden 1 , K. Otmianowska-Mazur 2 , H. Lesch 3 , and G. Skupniewicz 2 

1 Landessternwarte Heidelberg, Konigstuhl, D-69117 Heidelberg, Germany 

2 Astronomical Observatory, Jagiellonian University, ul. Orla 171, Krakow, Poland 

3 Universitats-Sternwarte, Scheinerstr. 1, D-81679 Miinchen, Germany 

Received; accepted 



ASTRONOMY 

AND 
ASTROPHYSICS 



Abstract. A fully three-dimensional computation of in- 
duction processes in a disk is performed to investigate the 
evolution of a large-scale magnetic field in gas-rich barred 
and spiral galaxies in the presence of field diffusion. As 
input parameters we use time dependent velocity fields 
obtained from self-consistent 3D N-body calculations of 
galactic dynamics. Our present work primarily analyzes 
the influence of the gas flows in non-axisymmetric grav- 
itational disturbances (like spiral arms and bars) on the 
behavior of the magnetic field. The magnetic field is found 
to evolve towards structures resembling dynamical spiral 
arms and bars, however its distribution is often much more 
complicated than the velocity field itself. We show that a 
random configuration of seed magnetic fields quickly dis- 
sipates, and is therefore unable to explain the magnetic 
field strength observed in nearby galaxies. The detailed 
comparison between the simulated velocity and magnetic 
fields allows us to establish that the reversals of magnetic 
vectors appear in the vicinity of the corotation radius of 
non-axisymmetric disk structures, where the velocity field 
shows rather abrupt changes. We argue that this explains 
the observed field reversals in our Galaxy. 

Key words: MHD - galaxies: magnetic fields - kinemat- 
ics and dynamics - spiral-galaxies 



1. Introduction 

The origin and evolution of large-scale magnetic fields 
in galaxies are still not properly understood (e.g. Beck 
et al. 1996; Lesch & Chiba 1997). There are several is- 
sues which are not yet clear. For example, there is still 
no consensus on the primordial origin of large-scale mag- 
netic fields as proposed by Kulsrud & Anderson (1992). 
The physical processes initiating and driving turbulent 
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galactic dynamos may not be able to excite large-scale 
fields because of the magnetic back reaction on the small 
scale turbulence (e.g. Cattaneo and Vainshtein 1991). Fur- 
thermore, the standard dynamo approach for investigat- 
ing large-scale magnetic fields does not include the non- 
axisymmetric part of the galactic velocity field which is 
able to explain several observed features in barred galax- 
ies (Chiba and Lesch 1994). We concentrate on this sec- 
ond issue in this contribution, which is an extension of 
our previous work (Otmianowska-Mazur et al. 1997, Paper 
I). In Paper I, we coupled, for the first time, fully three- 
dimensional velocity fields resulting from dynamical sim- 
ulations, with the large-scale magneto-hydro-dynamical 
evolution. Our main finding was an extremely sensitive 
response by the magnetic field to spatial variations of the 
gas velocity field. We concentrated on barred galaxies with 
normal gas content, where the magnetic vectors are found 
to show opposite directions along the sides of the bar, and 
in the course of disk rotation the field vectors are mixed 
at the ends of the bar. We have also found that the mag- 
netic field lines are well aligned with the bar and the spiral 
arms. 

In the present paper, we generalize our study to con- 
sider gas-rich spiral galaxies. The influence of the gas mass 
and gas contents on the structure and evolution of the 
magnetic field is discussed. 

In the magnetic field simulations, we are interested in 
galaxies which undergo significant dynamical evolution in 
terms of bar and spiral arm excitation. It is well known 
that the intensity of the corresponding dynamical insta- 
bilities depends sensitively on the amount of gas in the 
disk. The higher the gas content the faster and stronger 
are the excited non-axisymmetric gravitational instabili- 
ties which drive the spiral arms (Lynden-Bell & Kalnajs 
1972), as can be observed in actual galaxies (Rix & Zarit- 
sky 1995). It is also known that the excitation of such non- 
axisymmetric structures in disks leads to strong angular 
momentum transport toward the galactic center (Roberts 
et al 1979, Toomre 1977, Fuchs & von Linden 1997). 
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The angular momentum is either removed from or trans- 
ferred to the stellar population. This process develops the 
non-axisymmetric gravitational potential (Schwarz 1981). 
Since magnetic fields are strongly coupled with gaseous ve- 
locity fields, especially the velocity gradient, the amount 
of gas in the disk will influence the results for the mag- 
netic field evolution as well. In order to study this influ- 
ence in the present work, we perform calculations with 
gas masses of 10% and 20% of the total galactic mass. 
The simulation of a gas-rich galaxy allows us to focus also 
on proto-galaxies which have not undergone a complete 
transformation of primordial gas into stars. 

2. Models and numerical methods 

Massive galactic disks are able to form strong elongated 
bars (Sellwood 1981). This allows us to concentrate on 
a highly anisotropic velocity distribution accompanied by 
the development of a bar. Since the growth time scale of 
a bar increases with the ratio of the halo masses and the 
disk masses Mh/Md (Combes & Sanders 1981) we follow 
a long-lived strong bar by choosing Mh+b/Md = 2 (Md= 
stellar disk mass, Mh+b— halo and bulge mass). To fulfill 
the bar instability criterion 

% > (fl(r) - «(r)/2), (1) 

where Q p is the pattern velocity of the bar, £l(r) is the an- 
gular velocity of the disk and k{t) the epicyclic frequency. 
The disk has to be massive enough so that the mass in the 
inner part of the disk reaches a critical mass M > M cr i t 
(Sellwood 1981). We compare our simulations of a gas- 
rich barred galaxy (model III) with a model presenting 
only spiral activity in a less massive disk (model IV). 

2.1. N-body numerical scheme 

Our simulations are performed with a 3D N-body code in- 
volving a molecular cloud scheme as described by Combes 
& Gerin (1985) and Gerin et al. (1990). The modeled 
galaxy is embedded in a spherical halo consisting of visible 
and dark matter and in a bulge component in the central 
part of the disk. The galactic disk consists of gas clouds 
and stellar particles which are interacting via gravity. 

The stars interact only by softened gravitational forces, 
where the softening length e is much smaller than the scale 
height of the disk so that there is practically no artificial 
stabilization caused by large e parameters (Romeo 1994). 
The gas clouds can interact inelastically, which is simu- 
lated in an elaborate cloud-in-cell scheme describing the 
coalescence and fragmentation of the clouds. We adopt 
the same scheme and local processes as described by Ca- 
soli & Combes (1982) and Combes & Gerin (1985). The 
clouds encompass masses between 5 • 10 2 to 5 • 10 5 M 
distributed in 10 mass-bins with logarithmic intervals in 
order to take into account the power-law variation of the 
mass spectrum. 



Alternative methods to simulate the gas component as 
a fluid are smoothed particle hydrodynamics (e.g. Friedli 
& Benz 1993) or finite-difference schemes. Describing the 
dynamical features in disks, there is no qualitative differ- 
ence between these codes. In a future work, we will use 
different schemes for the gas particles to check the influ- 
ence on the magnetic field. 

We perform our simulation with N s ~ 38 000 stellar 
particles. To investigate the influence of the gas particle 
number on the galactic evolution we perform simulations 
with N g < 38000 and < 19000. We use the same grid size 
and FFT method for the calculation of the gravitational 
potential as described in Paper I. 

The time step St (=10 6 yr) of the N-body simulation is 
chosen to fulfill the relation v p hi, max 5t < I, where v p hi,max 
is the maximum of the rotation curve and / is the grid 
length (Combes & Sanders, 1981). The time step for the 
cloud-cloud collision is about 10 7 yr. Several different sim- 
ulations have been performed, varying the input param- 
eters for the disk, the bulge, the number of gas particles 
and the gas mass. 

2.2. Magnetic field evolution model 

We analyze the time dependent solutions of the induction 
equation: 

dB/dt = iot(v x B) - rotfarotB), (2) 

where n is the magnetic diffusion coefficient, v is the ve- 
locity of the gas and B is the magnetic field. As in Paper 
I, the computations are performed with the ZEUS3D code 
(Stone & Norman, 1992) using only that part which deals 
with magnetic field evolution. In the original form of the 
code the constrained transport (CT) algorithm is imple- 
mented for the flux-freezing approximation of the eq. (|2]) 
(7/ = 0). In our analysis, the effects of the magnetic diffu- 
sion are introduced in two ways: directly into the code us- 
ing the CT algorithm and by " pulsed flow" method (Bayly 
& Childress, 1989) described in Paper I. Both methods 
fulfill the condition divB = 0. 

To incorporate the particle-based velocity field into the 
grid-based code, a spline function interpolation is applied 
(cf. Paper I). 

The computations have been performed using Cray Y- 
MP in HLRZ Jiilich and Convex Exemplar SPP-1000 Se- 
ries in ACK CYFRO NET-KRAKOW. 

3. Input parameters 

We study numerically two experiments provided by our 3D 
particle-mesh scheme simulations: the barred galaxy (case 
III) and the spiral galaxy (case IV). For both the halo and 
the bulge components, Plummer spheres are adopted. In 
all models, the mass and the scale length of the dark mat- 
ter halo are the same as in Paper I. The model parameters 
for the N-body calculations are given in Table |. 
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Model 


III 


IV 


Mass in 1O 1O M : 






- dark-halo mass Mh 


9.6 


9.6 


- disk mass Add 


7.2 


5.8 


- gas mass M g 


1.6 (0.8) 


1.6 


- bulge mass Mb 


4.8 


6.2 


- mass ratio (M d + M g )/M to t 


0.38 (0.36) 


0.28 


Scale length in kpc: 






- dark halo h 


15 


15 


- bulge b 


1.1 


1.2 


- disk d 


7 


6 



Table 1. Input parameters for model III and IV 



Model 


dx, dy 


dz 


N x ■ 


Ny 


N z 


V 


Input 




[pc] 


[PC] 








[cm 2 /s] 


B-field 


Ilia 


300 


50 


101 


• 101 


■ 21 





toroidal 


Illb 


300 


50 


101 


■ 101 


• 21 


3.2 ■ 10 26 


toroidal 


IIIc 


300 


50 


101 


■ 101 


■ 21 





random 


Hid 


300 


50 


101 


• 101 


■ 21 


3.2 ■ 10 26 


random 


IVa 


300 


50 


101 


■ 101 


■ 21 





toroidal 


IVb 


300 


50 


101 


• 101 


• 21 


3.2 ■ 10 26 


toroidal 



Table 2. The magnetic field evolution model parameters 



Fig. 1. Distribution of molecular clouds for the model III 
projected onto the galactic plane (at t= 3.5- 10 s yr.). Mas- 
sive clouds are plotted with thicker dots as less massive 
clouds. 

The calculations for the magnetic field evolution are 
performed in 3D rectangular coordinates, where the XY 
plane is the galactic plane and the Z axis is the axis of 
galactic rotation. The rectangular size in the X and Y 
directions is 30 kpc. As a grid interval, we choose the dis- 
tance of 300 pc resulting in 101 grid points along the X 
and Y axis. The experiments performed in Paper I have 
revealed that simulations with twice the distance between 
the grid points are significantly affected by numerical dis- 
sipation due to truncation errors present in the ZEUS3D 
code (see Paper I, 3.2 and 4.3). Therefore, we decide to 
adopt a smaller interval which allows us to diminish the 
numerical diffusion to some extent. The scale height of our 
model galaxy is 1 kpc with 50 pc as the grid interval (21 
grid points along the Z axis). The model input parameters 
for the magnetic field evolution are summarized in Table 
||. The adopted time step of 10 5 yr fulfills the Courant 
condition for the given values of the velocity field and the 
magnetic diffusion coefficient. 

The computations were performed for two initial mag- 
netic field configurations: a purely axisymmetric toroidal 
field, going in the anti-clockwise direction, and a random 
field. The strength of the axisymmetric field is 1 /iG and 



-i 




-10 10 

[kpc] 



Fig. 2. Isophotes of the stellar disk density for the model 
III (at t= 3.5 ■ 10 s yr.). 

decreases both toward the galaxy center and toward the 
edges (see Paper I for a more detailed description). The 
random field is calculated numerically from a random dis- 
tribution of the magnetic vector potential A. Such a pro- 
cedure provides us with the magnetic field structure free 
of magnetic mono-poles. The mean strength of the ran- 
dom field is about 1/xG. As in Paper I, we adopt "out- 
flow" boundary conditions as the most stable for our cal- 
culations. The assumed turbulent diffusion coefficient is 
3.2 • 10 26 cm^" 1 in the XY plane and about 1.2 ■ 10 26 
cm 2 s~ 1 in the Z direction for the experiments Illb, Hid 
and IVb (see Table 2). The computations checking the in- 
fluence of different physical magnetic diffusion values on 
the magnetic field evolution were made in Paper I. The 
experiments III and IV last 6 • 10 8 yr. 

To compare the adopted numerical algorithms for the 
magnetic diffusion term in the induction equation (j^) an 
initially random field is evolved without any velocity field 
(« = 0). 

4. Pattern velocity and Lindblad resonances 

The non-axisymmetric structure in the disk has an angu- 
lar pattern speed given by the angular velocity Q p . Three 
radii in each galaxy are of specific interest because the 
pattern speed resonates with the eigenfrequencies of the 
unperturbed galaxy. The first order epicyclic theory (e.g. 
Binney and Tremaine 1987) superimposes harmonic os- 
cillations on the rotation in both the radial and tangen- 
tial directions with a characteristic frequency called the 

1 /2 

epicyclic frequency n = 2ft [1+ . Thus, in the 

proper rotating frame, the particle will move in a ret- 
rograde sense around a small ellipse, called an epicycle, 
with axial ratio K/2ft. The resulting motion in the in- 
ertia! frame is a rosette orbit, generally not closed. To 
the first approximation, the motion of a gas cloud or a 
star is thus governed by two frequencies. In a frame coro- 
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Fig. 3. Number density of clouds per mass bin (logarith- 
mic mass intervales) for the model III at t = 1.5 x 10 8 yr 
(solid line) and t == 6.0 x 10 8 yr (dotted line). Mass bin 1 
is the lowest cloud mass class ( 5- 10 2 M©); mass bin 10 is 
the most massive one (5- 10 5 M©) 

fating with a spiral these are fl(r) — Cl p (r) and n(r). 
When they are comparable, i.e. when the relative fre- 
quency v == (u> — mO)/« == m(Ct — Q, p )/k is equal to a 
rational number, we have a resonance. The three main res- 
onances in a galaxy (with a two-arm spiral pattern m = 2) 
are 

- the inner Lindblad resonance (ILR) where v = — 1 
and f2(r) — n(r)/2 == Q, p 

- the corotation resonance (CR) where v = and 
n(r) == n p (r) 

- the outer Lindblad resonance (OLR) where v = +1 
and fi(r) + n(r)/2 == tt p 

At the Lindblad resonances, the disk and the pat- 
tern (bar, spiral arms) intensively exchange energy, mass 
and angular momentum. The radial velocity, for exam- 
ple, changes its sign at CR. Inside CR, the material flows 
inward and outside CR it moves outward. Due to the en- 
ergy and angular momentum transfer, the orbits of the 
particles change by 90 degrees at every resonance. As we 
shall see, this simple picture is complicated by the fact 
that a bar as well as a spiral system has its own reso- 
nances. Due to the interaction of the particles with the 
non-axisymmetric structure, the bar and spiral arms are 
slowed down, thereby moving the resonance radii outward. 
Larger portions of the disk are then involved in the radial 
inflow of the material. Since the resulting disk velocity 
field, which we use as an input parameter for our mag- 
netic field calculation, reflects this dynamical behaviour, 
we analyze the pattern velocity and the radii of the Lind- 
blad resonances and corotation in our simulated galaxies. 

Based on a Fourier analysis of the density distribution 
or of the gravitational potential, which allows us to decom- 
pose the perturbation present in the galaxy (Junqueira & 
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Fig. 4. The rotation velocity and the radial velocity at 
t = 1.5 x 10 8 yr (solid line) and t == 6.0 x 10 8 yr (dashed 
line) for the model III 
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Fig. 5. The radial velocity at a 1 kpc thick disk stripe 
(z == 0) at different time step are shown. 

Combes 1996), we are able to calculate the pattern veloc- 
ity of our simulation at different time steps. By plotting 
the pattern velocity fl p together with the curves fl ± k/2 
we determine the radii of the resonances in the disk. 

The uncertainty of the pattern velocity amounts to 
about ± 2km/s/kpc. For the Lindblad radii the error in- 
creases from the inner to the outer part of the disk. 

5. Results of Model III 

5.1. 3D N-body simulation 

Simulations of barred galaxies were performed by several 
authors (see cf. Katz & Gunn 1991; Gerin et al 1990; 
Combes & Elmegreen 1993; Norman et al 1996). As in the 
present calculations, a strong spiral structure is formed 
first. The two-armed spiral configuration becomes more 
pronounced and more massive during the next ten time 
steps. After 1.5 • 10 8 yr this structure is clearly visible in 
the gaseous and stellar disks. 
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Fig. 6. a Magnetic field vectors and the gas density (grey plot) are plotted on the left for the model Ilia at t— 
1.5 • 10 8 yr. The belonging intensity of the magnetic field (grey plot) and the radii for the Lindblad resonances for the 
bar (CR^gj. and OLR^^., solid line) and for the spiral system ( CRsp, dashed line) are shown in the right figure (the 
ILR S p are at the same radius as the OLR^^; the OLR S p are located outside the plot range). 

Fig. 6. b As in Fig. ||a, at t = 2.5 • 10 s yr. The radii of the Lindblad resonances moved outwards during the simulation. 
Fig. 6. c As in Fig. ^|a at t = 3.5 • 10 8 yr. In the inner part of the disk a ILR^j. appears. 



In comparison with model I (see Paper I) , experiment 
III exhibits a higher mass concentration in the central part 
of the disk due to the assumption of a bulge of smaller scale 
length. This means that in this model, the bar appears 
earlier than in case I (after 2.3 • 10 s yr). A second reason 
for earlier bar formation is the larger amount of gas in the 
disk (see also Friedli & Benz 1993, Little & Carlberg 1991). 
After 2.8 • 10 s yr of evolution the bar gets thicker, and 
during subsequnt time-steps the gas concentrates along 
its major axis. A massive elongated bar is formed which 
ends near the corotation radius, as in the simulations of 
Sellwood (1981) and Combes & Sanders (1981). 

After 3.5 • 10 8 yr the corotation radius of the bar is 
located at 7.2 kpc. The bar has a length of 14.2 kpc and it 
concentrates 27% of the disk gas mass (further values are 
summarized in Tab|^). In the outer part of the disk, a two- 
armed spiral pattern is still visible (Fig.[l]). The pattern 
velocity and the Lindblad radii are given in Table § 

The orbits of the bar in the xy plane are similar to 
those of the 2D-simulation of Sparke & Sellwood (1987) 
and the 3D simulation of Hasan et al. (1993). The high 
central gas concentration leads to the appearance of X2 or- 
bits inside the ILR. Due to the thickness of the bar, these 
orbits have a large population. We show the isophotes of 
the stellar disk density in Fig.|[ It is clear that these X2 
orbits create the isophote twist visible in the inner part of 
the galaxy. The isophote twist is also present in the simu- 
lations of Friedly & Martinet (1993) (see e.g. NGC 1422) 
and it cannot be interpreted as a bar in bar (Shaw et al. 
1993). 

The mass spectrum of the molecular clouds changes 
during the evolution of the disk, which is plotted in Fig. || 
within 30 kpc of the disk. After 6 • 10 8 yr there are fewer 
particles in the disk; in particular the number of low 
mass clouds decreases because of cloud-cloud collision. 
The mass spectrum of the clouds in the bar is similar 
to the mass spectrum of the total disk. However, the per- 
centage of extremely massive clouds in the bar is smaller 
than that in the disk. 

During subsequent evolution, the pattern velocity of 
the bar slowly decreases and the bar shrinks as in other 
simulations, cf. Little & Carlberg (1991). 



5.2. The velocity field 

The velocity field of the disk given at different time steps 
is the input for magnetic field calculation. Fig. [| presents 
the rotation curve and the average radial velocity at dif- 
ferent time steps calculated in rings of constant distance. 
Negative v r vectors are directed toward the galactic cen- 
ter. Obviously, there is mass transport toward the center, 
as might be expected. Inside the bar, the radial veloc- 
ity is dominated by the £2-orbit family (Contopoulous & 
Papayannopoulos 1980) which allows for positive radial 
velocities. To examine the local radial velocity in more 
detail, a cut in the disk at z = through the midpoint 
with a width of 1 kpc is made for two time steps (Fig. ||) . 
At 6 • 10 s yr the cut is taken along the bar. This figure 
shows that the radial velocity reflects the dynamical be- 
haviour of the disk: e.g. at t=6-10 8 yr the ILRf, ar is located 
at 3 kpc and v r changes sign at this radius. 

5.3. The 3D evolution of magnetic field: model III 

Four simulations of barred gas-rich galaxies have been 
performed, charactaized by the inital configuration of the 
magnetic field: a toroidal magnetic field with (Illb) and 
without magnetic diffusion (Ilia) and a random field with 
(Hid) and without diffusion (IIIc) (Table ||). In Figs.|a- 
d the evolution of the large-scale magnetic field for the 
case Illb is presented. We show only the case with non- 
zero turbulent magnetic diffusion r\. However, the com- 
putations with rj = result in similar structures (but 
with a higher magnetic energy density, see Fig. ||) . For the 
present calculations, the diffusion coefficient rj is about 
three times smaller than in the case Id presented in Pa- 
per I, where diffusion significantly influenced the magnetic 
field behaviour. The models in Paper I are also affected 
by numerical diffusion. Thus, the present simulations su- 
persede those in Paper I. 

At the first evolutionary stage, the magnetic field 
quickly responds to the velocity shear of the interstel- 
lar gas flow in the non-axisymmetric gravitational distur- 
bances. Fig.^ja presents the magnetic field distribution at 
t = 1.5- 10 8 yr. The field vectors (left), shown in the galac- 
tic plane, are superimposed onto the gas density grey plot 
(in logarithmic scale). The magnetic field intensity distri- 
bution (right) is plotted together with the circles of the 
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Fig. 6. d As in Fig. ||a at t = 6.0 • 10 8 yr. The resonances for the spiral system are outside the plot range. 



timestep Q bar IL R D ar CR bar OLR bar 



[10 8 yr] 


r km I 
L skpc J 


fisp 

r km 1 
L s kpc J 


[kpc] 


[kpc] 


= ILRsp 
[kpc] 


CRsp 
[kpc] 


OLRsp 

[kpc] 


1.5 


66.9 


21.9 




4 


6.5 


13 


18.9 


2.5 


57.3 


14.3 




5 


8.1 


18.0 


23 


3.5 


31.8 


12.2 


2.1 


7.2 


11.5 


24.1 


27.3 


6.0 


25.0 


11.0 


3 


9.9 


17.9 


29 


~ 35 


± 


2 


2 


2.4 


2.4 


2.6 


3 


4 



Table 3. Pattern velocity of the bar and the spiral system (sp) and the radius of the corresponding resonances 



Linblad resonances of the the bar (CR^^ and OLR^^ 
solid line) and of the spiral arms (CRsp dashed line). 
From the beginning, the magnetic field closely follows the 
dynamical spirals resulting in a four-arm pattern visible 
in the figure. In one of the arms the orientation of the 
magnetic vectors is opposite to the initial anti-clockwise 
direction of those in the toroidal magnetic field. This is 
a result of the strong velocity shears present in this re- 
gion (see discussion in Sec. 7). There are also two regions 
with vector reversals located in the inner parts of the next 
two arms, close to the CR of the bar (see Figj^a). The 
magnetic intensity distribution (Fig.||a, right) shows that 
the magnetic features form nonuniform stripes exhibiting 
generally spiral structure. FigJ^b presents the magnetic 
field configuration at 2.5 • 10 8 yr. The circles of the reso- 
nances for the bar (CR ^, ar and OLR^^,, solid lines) and 
for the spiral arms (CRsp, dashed line) are superimposed 
onto the intensity grey plot (right). The dynamical struc- 
ture shows the initial phase of bar formation (see Fig.|^b) 
and two spiral arms. The magnetic field is aligned along 
the bar and spirals, but its distribution is more compli- 
cated. In the central part of the right arm (Fig.^Jb, left), 
the magnetic vectors are directed opposite to the vectors 
at the arm edges (Fig.^b, left), which is also visible in 
the intensity grey plot (Figj^b, right). In the upper part 
of the figure, the region with magnetic vector reversal is 
still present but the field is weaker than in those mag- 
netic features with an anti-clockwise direction of the vec- 
tors (Fig.^jb, right). In Figjfjc, we show the magnetic field 
distribution at t = 3.5 • 10 yr. The ILR now appears in 
the bar. Three resonances (ILRj-^, CRj-^j. and OLR^^) 
for the bar are marked here by solid lines. In the outer 
part of the disk, the spiral arm system is still visible (see 
Fig.^). The magnetic field vectors are again aligned with 
these structures but also form additional arms and regions 
with vector reversals. As in experiments I and II of Pa- 
per I, the magnetic field vectors start to be mixed at both 
ends of the bar. In FigJ^c, magnetic spiral features form 
a pattern of four non-uniform spirals (see Fig^Jc, right). 
However, there is now a considerable difference between 



Fig. 7. Topolines of the radial velocity superimposed onto 
the magnetic intensity grey plot for the simulations Illb 
at t = 1.5 ■ 10 8 yr 

the present and 1.5 ■ 10 8 yr time step (Fig.||a). In the pre- 
vious evolutionary stage, magnetic spiral features closely 
follow the dynamical spirals. At t = 3.5 ■ 10 8 yr, the mag- 
netic field has two more arms than the gas (Fig. [3)3). The 
regions with vector reversals are still present between the 
two left magnetic spiral features, along the right side of 
the bar and in the inner part of the right arm (Fig.||c;). A 
broad bridge of magnetic vectors connects the lower end 
of the bar with the right side arm and it has no significant 
dynamical counterpart. 

The main conclusion from these facts is that the 
magnetic field "remembers" the dynamical dis- 
turbances significantly longer than they exist. At 
t = 6 • 10 8 yr, the situation is even more evident. The 
dynamical arms disappear but magnetic spirals remain 
(Fig.||d). The magnetic field lines are aligned with the 
bar and create two diffused lanes consisting again of two 
sub-arms with reversed magnetic vectors (Fig.|6jd, left). In 
the right arm, the magnetic intensity is similar to stripes, 
both, with and without reversal (Fig.||d, right). 

The evolution of the galactic magnetic field under the 
influence of non-axisymmetric gas flows shows that the 
magnetic field reacts quickly even to weak disturbances in 
the field. The direction of the magnetic vectors can eas- 
ily be changed by velocity shear, especially in the regions 
close to the Lindblad resonances of the bar and in the spi- 
ral arms. We could connect these field reversal with the 
corotation radius of the disk (see discussion in Sec. 7). The 
numerical simulations also show that magnetic field vec- 
tors are mixed at the bar ends, as was shown in Paper 
I. 

In order to analyze the dependence of the magnetic 
field intensity distribution on the velocity field shear, FigJ^ 
presents the topolines of the radial velocity v r over the 
magnetic intensity grey plot for the simulations Illb at 
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Fig. 8. Time evolution of the magnetic energy density E 
normalised to the initial energy density E Q for the models 
III. The figure legend is in agreement with Table 2 

t = 1.5 • 10 8 yr. The maximum of the magnetic intensity 
is located in the region with the strongest gradient of the 
v r component, which coincides with the CR^gj. circle (see 
Fig.^a). On the inside of the main magnetic spiral fea- 
tures (upper-central part of the figure) the velocity distri- 
bution has a maximum, while, on the outside, a minimum 
is visible (see Fig|j]). Both maxima are distributed along 
the OLR^gj resonance circle (Fig.||a) . The next magnetic 
arm, which is seen in the lower part of the figure, is also 
connected with such reversals (CRsp radius) of radial ve- 
locity. The magnetic strength of this arm is weaker than 
that of the previous one, but the maxima and minima on 
both sides of the arm are not so extreme. This supports 
the conclusion that the magnetic field reversals ap- 
pear due to the reversal of the velocity at the CR 
radii. The rest of the magnetic features are not so promi- 
nent and it is difficult to associate them with changes in 
the velocity shear. 

Fig. |^ presents the time evolution of the magnetic en- 
ergy density (E) normalised to the initial energy density 
(E ) for the four cases in model III. The parameters if 
the four vases are given in Table 2. The energy density 
curves show the time evolution of the toroidal (case Ilia 
and b) and random (case IIIc and d) input magnetic field 
structure. Simulation Ilia, with no physical diffusion (solid 
line) , shows a constantly decreasing magnetic energy den- 
sity due to numerical effects and the lack of strong ve- 
locity shears compared to the barred galaxy case in Pa- 
per I. However, three small maxima are still present (due 
to the spiral arms and the bar formation). In model Illb 
(dashed-dotted line) the turbulent diffusion diminishes the 
magnetic energy density much faster than in the previous 
case. The density for simulation IIIc (short dashed line) 
with the random input field configuration increases by a 
factor of 2.5 and then preserves the energy density value 



during the whole evolution. This maybe explained by the 
random character of the initial magnetic field and/or the 
lack of physical diffusion. The characteristic scale of mag- 
netic field variation has a length of one grid step, so the 
numerical resistivity is not efficient enough to diminish the 
field. 

The problem of the numerical resistivity was tested in 
several experiments with different input random field and 
without physical diffusion. The magnetic energy decreased 
when the input random magnetic field possessed the scale 
of change larger than one grid step. The random input 
magnetic field with magnetic vector randomly distributed 
at each grid point avoid numerical resistivity. This fact is 
probably connected with the numerical resistivity which 
is not strictly a magnetic one (Hawley & Stone, 1995). 
Stationary fields, for example, do not diffuse away. 

All other experiments made with the random magnetic 
field and with the turbulent diffusion exhibit a rapid de- 
crease in the magnetic energy density which is present in 
case Hid (long dashed line). The random magnetic vectors 
are distributed in all directions and, under the influence 
of the diffusion, they easily interact, diminishing the mag- 
netic field strength. The decrease of magnetic energy den- 
sity in each experiment strongly suggests that, in order to 
explain observed galactic fields, a physical mechanism (e.g. 
a magnetohydrodynamical dynamo) is necessary for mag- 
netic field amplification on the scale of the whole galaxy. 

In order to compare the numerical schemes used in our 
experiments (the diffusion is implemented directly into the 
code with the CT algorithm and with the " pulsed-flow" 
method (see Paper I)), we perform two experiments with 
random input magnetic field structure and with no veloc- 
ity field. Both methods give a similar slope of decreas- 
ing magnetic energy density. However, the CT scheme 
maintains a slightly higher magnetic energy density than 
the "pulsed-flow" method. This is because the "outflow" 
boundary conditions used in the calculations with the CT 
method, which are otherwise the most stable for our sim- 
ulations, exhibit small magnetic monopole values at the 
boundaries. 

6. Results of Model IV 

6.1. 3D N-body simulation for model IV 

Model IV also assumes a developing two-armed spiral pat- 
tern. This spiral structure transports angular momentum 
into the outer part of the disk and, outside the corota- 
tion radius CRsp, the gas particles gain angular momen- 
tum. Due to the assumption of a less massive disk in this 
model, no bar appears in the disk during the simulation. 
The mass concentration in the inner part of the galaxy 
is maintained during the evolution, but the arms become 
weaker and weaker (after approximately 2.5 • 10 s yr). A 
spiral structure is sharpest at the earliest stages of evo- 
lution (see also the simulation in Elmegreen & Tomasson 
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Fig. 9. a Magnetic field vectors and the gas density (grey plot) are plotted left for the model IVb at t= 1.0 • 10 8 yr. 
The belonging intensity of the magnetic field (grey plot) and the radii for the Lindblad resonances for the spiral system 
(solid line) are shown in the right figure. 

Fig. 9. b As in Fig.fJa now for t=1.5 • 10 s yr. The positions of Linblad resonances move outward. 



timestep Qsp ILRsp CRsp OLRsp 



[10 s yr] 


r km I 
L skpc J 


[kpc] 


[kpc] 


[kpc] 


1.0 


28.6.9 


2.5 


9.0 


15 


1.5 


24.31 


3 


9.8 


17.8 


2.0 


22.16 


3.5 


11.3 


19.3 


2.5 


18.59 


4.2 


13.5 


26 


± 


2 


4 


4 


4 



Table 4. Pattern velocity the spiral system and the radius 
of the resonances for model IV 



1993). The vertical scale of the stellar disk changes dur- 
ing evolution from ± 0.5 kpc (1 • 10 8 yr) to ± 1.5 kpc 
(2.5 • 10 s yr), especially at the outer parts of the disk. The 
pattern velocity and the resonances of the spiral mode are 
calculated as described above. At 1.5 ■ 10 8 yr, we find the 
pattern velocity of the m=2 mode is about 25 km/s/kpc 
(see Table ^) . The spiral structure is visible for the next 
4 • 10 s yr. 

6.2. The 3D evolution of magnetic field: model IV 

The evolution of the magnetic field in the spiral galaxy 
(model IV) is caused by an efficient and fast expulsion of 
the magnetic field from the central part of the galaxy and 
is followed for 4.0 ■ 10 s yr. 

The magnetic field vectors superimposed onto the mass 
density distribution (grey plot, left) and the magnetic field 
intensity distribution (grey plot, right) shown together 
with the Lindblad radii of the m=2 spiral mode system 
(solid circles) are presented in Fig.^a-c at four times: 1.0, 
1.5, 2.0 and 2.5 ■ 10 8 yr for the case IVb. 

At 1 ■ 10 8 yr (Fig.^a), the magnetic field forms two dif- 
fuse spiral arms. In the arms, the magnetic intensity is 
non-homogeneous. The areas of stronger and weaker field 
density are visible, though outside the arms the intensity 
distribution is more uniform. The field is quickly expelled 
out of the central part of the galaxy, which possesses a 
strong mass concentration. This is the case throughout 
the whole evolution for 4.0 ■ 10 8 yr. At this early stage, 
the bifurcations of the magnetic spiral features are ob- 
served. The magnetic field vectors are aligned along the 
outer parts of the visible arms and regions of magnetic 
field reversal become visible (FigJ^a, left). For 1.5 • 10 s yr, 
the magnetic spiral features become more distinct and 
characteristic " islands" of stronger fields are visible along 



them. Thus, the field is still inhomogeneous and is driven 
out from the galactic center. It forms a spiral pattern. The 
magnetic field vectors are aligned along visible arms and 
reversals are stronger and more clearly visible. Areas with 
vectors in opposition to the initial toroidal field are again 
found at the CR circle (see Fig.^ja, right). For 2.0 • 10 8 yr, 
the field is driven out of the central part of the galaxy, 
and the magnetic spiral features are nonuniform. They 
are clear and sharp at the beginning and more diffused 
at their ends. They begin to merge with each other so 
that the reversals almost disappear. The magnetic spirals 
extend further out of the disk than their dynamical coun- 
terparts. The magnetic field intensity starts to develop 
a ring at the corotation radius. At 2.5 • 10 8 yr, the spi- 
ral structure of the magnetic arms disappears. The field 
forms a clear ring outside the dynamical structure, still 
aligned with the CR circle (see discussion in Sec. 7). Its 
vectors are in the anticlockwise direction and possess an 
almost azimuthal component (axisymmetric spiral struc- 
ture - ASS?). It is still inhomogenous with "islands" of 
stronger field density. Very weak magnetic spiral features 
appear in the center, probably caused by the numerical 
artifacts. 

The magnetic energy density curves for model IV are 
presented in Fig.|l^. The experiment with no physical dif- 
fusion (IVa, solid line) results in two maxima of magnetic 
energy density near 0.7 • 10 s yr and 1.1 • 10 8 yr. The first 
maximum is connected with a strong velocity shear due 
to the initial formation of spiral arms in the galactic disk. 
The second maximum is the result of the velocity gradient 
caused by a strong mass transport into the central part of 
the galaxy. The galaxy develops a massive spiral structure 
that extends to the centre of the disk. After two maxima, 
the energy density decreases due to numerical resistivity, 
but still possesses further maxima caused by continuous 
changes in the spiral pattern (about 1.6 • 10 s yr) and later 
at about 2.9- 10 s yr and 3.6- 10 s yr. These two maxima are 
connected with the formation of a small bar in the central 
part of the galaxy. The magnetic energy density for the 
case IVb with physical diffusion (dashed line) decreases 
for the whole evolution showing only small maxima at the 
same positions as in previous simulations. 

The analysis of B z , topoline maps made at different 
heights of the galactic disk, shows that weak magnetic 
field waving is present in both models (III and IV) (see 
Paper I for discussion). 
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Fig. 9. c As in Fig.|a now for t=2.0 • 10 8 yr. 

Fig. 9. d As in Fig.|]a now for t=2.5 • 10 8 yr. The ring for the OLR is now outside the plot range. 
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Fig. 10. Time evolution of the magnetic energy density E 
normalised to the initial energy density E Q for the model 
IV 

7. Conclusions and discussion 

A numerical model involving the three-dimensional and 
time-dependent velocity of molecular gas is applied to 
study the evolution of the large-scale magnetic field under 
the influence of non-axisymmetric dynamical structures. 
We perform a comparison between bar and spiral arm ac- 
tivity to investigate the effect of the different gas velocity 
fields. The bar transports much more mass inward and an- 
gular momentum outward than the disk with spiral arms 
alone. Our present work demonstrates that the magnetic 
field is influenced more strongly by the bar velocity shear 
than by spiral disturbances. Our main conclusion is: 

• The magnetic field dissipates much faster in spiral than 
in barred galaxies due to the lower level of angular 
momentum transport. 

The number of particles in the gas disk has no significant 
influence on magnetic field behaviour. 

In all experiments, the regions with reversals of the 
magnetic vectors appear due to velocity shears. We note 
that: 

• Field reversals appear close to the corotation radius 

at which, according to standard theory (e.g. Athanassoula 
1984), the radial velocity of the gas flow changes its sign. 
Outside the CR circle, the radial velocity is directed out- 
ward while inside the CR, the gas flows inward if the disk 
has only one resonance system. At the CR, there is a 



radial-velocity minimum. Since the magnetic field is cou- 
pled to the gradients of the velocity field, it changes its di- 
rection as well. In addition to the corotation radius of the 
bars, there is also a corotation radius of the spiral arms. 
Similarly, at the CRsp, the magnetic field changes its di- 
rection in response to the varying velocity field at that 
radius. Such behaviour can be observed in all our simula- 
tions but especially in model IV (spiral galaxy) where the 
obtained ring of magnetic field (Fig. 9c & d) shifts together 
with corotation radius of the spiral arms. 

Changes in the velocity field at the inner Lindblad res- 
onance (ILR) are also visible in the magnetic field. The 
orbits change their sign by 90 degrees at every resonance. 
Superimposed onto the global galactic rotation, this leads 
to a large-scale spiral pattern, whose complicated veloc- 
ity structure is traced by the magnetic field via velocity 
gradients. 

From our numerical experiments, we demonstrate that 
due to a smaller number of Linblad resonances in spirals 

• the regions of magnetic vector reversals are less signif- 
icant in spiral than in barred galaxies. 

In addition, we have studied how the input magnetic 
field affects the further development of the magnetic field 
distribution and its intensity. We have examined purely 
axisymmetric toroidal fields with an anti-clockwise direc- 
tion and random seed fields. 

The case of non-axisymmetric input magnetic field has 
been widely discussed in the paper Otmianowska & Chiba 
(1995). We analyzed the uniform protogalactic magnetic 
field, distributed along the Y axis. For all experyments, 
the magnetic field was quickly wound into the bisymmet- 
ric structure, causing the initial growth of the magnetic 
energy. As the time progresses the magnetic field in the 
inner part of the disk was quickly swept out due to mag- 
netic lines winding mechanism and the magnetic energy 
always decreased very fast (after 2 x 10 s yr). 

In contrast to the computations with a toroidal mag- 
netic field, those with random magnetic fields resulted in 
an extremely rapid decrease in the magnetic energy den- 
sity. This suggests, that: 

• Random seed magnetic fields in galaxies cannot ex- 
plain the observed magnetic field in galactic disks. 

As an interesting application, we note that observa- 
tions of the magnetic field in our Galaxy show strong ev- 
idence for two reversals of the field vectors (Vallee 1988, 
1991, Han & Qiao 1994). Vallee (1997) argues that the 
region with magnetic vectors directed opposite to the rest 
magnetic field in our Galaxy (which has a clockwise direc- 
tion) extends from the Sagittarius arm (about r — 6.5 kpc) 
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to the location of the Sun (7.5 ± 1 kpc). These observed 
reversals of magnetic field are located in the vicin- 
ity of the observed corotation radius in our Galaxy. 

Mishurov et al. (1997) estimated the pattern velocity of 
the spiral pattern as £1 = 28.1 ±2.0, and this sets the coro- 
tation radius of our Milky Way at a radius of 7.2 ±1.3 kpc. 
This might confirm the statement that the magnetic field 
vectors' distribution and direction depend on the dynam- 
ical behaviour of the non-axisymmentric structures and 
the positions of the Lindblad resonances as well. Existing 
theoretical explanations (e.g. Poezd et al. 1993) concern- 
ing the observed reversals in our Galaxy connect them 
with special conditions of thin dynamo action and with 
structures inherited from seed magnetic fields. Our results 
provide a simpler solution to this problem. The regions 
of magnetic reversal appear due to clear dynamical pro- 
cesses, closely connected with non-axisymmetric features 
present in a galactic disk. In a future work we will in- 
vestigate whether there are more such "coincidences" in 
other galaxies. For example, the ring structure in model 
IV might be interesting for our neighbouring galaxy M31 
in which such a ring was detected 15 years ago by Beck 
(1982). 

However the modeled ring is shifted further from the 
galactic center than the gaseous arms visible in Fig.9d. 
The observed magnetic ring coincides with the HI torus 
at the radius of 10 kpc. 

The field reversals are also interesting for comparison 
with the observations of M81 by Krause et al. (1987a, b), 
who determined a neutral line at which the magnetic field 
strength should vanish. 
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